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In this paper we revisit and update the computation of thermal corrections to the stability of the 
electroweak vacuum in the Standard Model. At zero temperature, we make use of the full two-loop 
effective potential, improved by three-loop beta functions with two-loop matching conditions. At 
finite temperature, we include one-loop thermal corrections together with resummation of daisy 
diagrams. We solve numerically—both at zero and finite temperature—the bounce equation, thus 
providing an accurate description of the thermal tunneling. Assuming a maximum temperature in 
the early Universe of the order of GeV, we find that the instability bound excludes values of the 
top mass Mt > 173.6 GeV, with Mh — 125 GeV and including uncertainties on the strong coupling. 

We discuss the validity and temperature-dependence of this bound in the early Universe, with a 
special focus on the reheating phase after inflation. 


I. INTRODUCTION 

The Standard Model (SM), if extrapolated up to ex¬ 
tremely high energies by means of its Renormalization 
Group (RG) equations, reveals a rather peculiar prop¬ 
erty: the electroweak vacuum does not correspond to 
the configuration of minimal energy; contrarily, it is a 
metastable state close to a phase transition [I]. This 
scenario—dubbed near-criticality —must be considered 
as the most important theoretical message learned from 
the LHC run 1. Near-criticality may open a window on 
the realm of Planck-scale physics, otherwise completely 
inaccessible from a phenomenological point of view. Un¬ 
derstanding its meaning, and refining the computational 
tools needed to this end, is therefore a task of primary 
importance. 

On the quantitative level, near-criticality emerges from 
the computation of the tunneling probability—integrated 
over the age of the Universe—between the false and 
true vacuum of the Higgs potential [2, 3]. Probabili¬ 
ties larger than one correspond to an unstable configura¬ 
tion of the electroweak vacuum. In the SM with Higgs 
mass Mh — 125 GeV the instability occurs if Mt > 178 
GeV, a value of the top mass that is fairly away from 
present experimental measurements. In other words, for 
the present central values of Mh and Mt , the electroweak 
vacuum of the SM is unstable but sufficiently long-lived 
if compared to the age of the Universe. 

However, this result relies on the assumption that ther¬ 
mal effects, due to non-zero values of the temperature, 
are neglected. The impact of thermal corrections on the 
computation of the tunneling probability was intensively 
discussed in the past [4-6]. Intuitively, thermal fluctua¬ 
tions at finite temperature increase the tunneling prob¬ 
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ability, and the easiest way to visualize their role is to 
think about the analogy with the one-dimensional quan¬ 
tum mechanical system of a particle in a potential with a 
false ground state. The thermal kinetic energy borrowed 
from the heat bath shifts the particle from the initial po¬ 
sition at the bottom of the false vacuum, thus facilitating 
the tunneling across the potential barrier. In quantum 
field theory a proper formulation of the problem requires 
the computation oi i) the finite temperature effective po¬ 
tential, and ii) the bounce field conhguration, namely the 
solution of the classical equations of motion that triggers 
the tunneling between the false vacuum and the other 
side of the potential barrier [2, 3]. 

Apart from computational technicalities, on the inter- 
pretational side thermal corrections may play an impor¬ 
tant role since it is very likely that our Universe—in the 
early stages of its existence—went through an extremely 
hot phase. In [7, 8] the instability of the electroweak 
vacuum was investigated from a cosmological perspec¬ 
tive (see also [9-16]). The main emphasis of [8] was put 
on the computation of quantum fluctuations of the Higgs 
field during inflation. The heart of the matter is that 
these fluctuations may force the Higgs field to fall down 
into the true minimum even before inflation ends. How¬ 
ever, the bottom line is that this problematic situation is 
not realized if the reheating temperature after inflation 
is sufficiently large. As a consequence, [8] points towards 
a cosmological scenario in which, right after inflation, 
the Universe is characterized by an extremely high value 
of the temperature. Under this condition, thermal cor¬ 
rections to the tunneling probability—as already noticed 
in [8]— can not be neglected. 

Motivated by this result, in this paper we revisit and 
update the computation of thermal corrections to the 
stability of the electroweak vacuum in the SM, and we 
structure our work as follows. In section H, we discuss the 
finite temperature effective potential used in our analy¬ 
sis. In section HI, closely following the approach of [6], 
we compute the bounce solution and the probability of 
thermal tunneling. In section IV, we present our results 
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in terms of the so-called phase diagram of the SM. Fi¬ 
nally, we conclude in section V. In appendix A, we pro¬ 
vide further detail about the finite temperature effective 
potential presented in section II. 


II. EFFECTIVE POTENTIAL AT FINITE 
TEMPERATURE 

As stated in the introduction, the starting point of 
our analysis is the effective potential of the SM at finite 
temperature. We use the following short-hand notation 

Ueff((^,T) = UoW+Vi_loop(</>) + 'F2-loop(</>) 

+ Vi-loop(</>,r) -h Vcing (1) 

where the first (second) line refers to T = 0 (T ^ 0), 
and (j) is the real Higgs field. At T = 0, we include, 
in addition to the tree level Higgs potential Vo{(j)), one- 
and two-loop corrections. At T 7 ^ 0, we include one- 
loop thermal diagrams and plasma effects, the latter de¬ 
scribed by one-loop ring resummation of daisy diagrams. 
For completeness, we collect the explicit expressions in 
appendix A. In appendix B we discuss the validity of 
the one-loop approximation at finite temperature. On a 
more technical level, the effective potential in eq. ( 1 ) is 
equipped with the following tools. 

o We implement the RG improvement of the effective 
potential in eq. (1). The dimensionless parameters 
run according to the three-loop SM RG equations. 
The running Higgs field is 

(j){t) = , r(t) = f , (2) 

Jo 

with 7 (t) the Higgs field anomalous dimension 
d(j){t)/dt = 'y(t)cj){t). 

o The matching condition are evaluated at two loops, 
following [I]. 

o In order to canonically normalize the Higgs ki¬ 
netic term, we introduce the canonical field 0 can = 


^ After renormalization the Higgs field effective lagrangian is 

tleff = , (3) 

with r(^) = The canonically normalized Higgs 

field 0can is implicitly defined by d^can/d^ = We use 

the approximate solution (jjcan — This approximation 

amounts to take a constant It corresponds to 7 ( 0 ) 1, 

since d(f)canld(p = + 7(0)]. Indeed, we checked that this 

condition is always verified during the RG evolution. 


o Finally, in order to minimize the impact of large 
logs, the renormalization scale is chosen according 
to 

= 7an , (4) 

where the relation with the running parameter t is 
fi{t) = /ioexp(t). The scale /xq fixes the starting 
point of the running, and we use as a reference the 
physical top mass. From now on, we suppress the 

subscript can- 


M, = 173.34 GeV, M* = 125.09 GeV, as = 0.1184 
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FIG. 1: Effective potential and first derivative at zero and 
finite temperature as a function of the Higgs field. Blue line: 
T = 0. Red line: finite temperature T = lO'^® GeV. Solid line: 
effective potential. Dashed line: first derivative. The values 
of the input SM parameters are shown in the plot label. 

We show the effective potential at zero (blue solid 
line) and finite (red solid line) temperature in fig. 1. 
Dashed lines refer to the corresponding first derivative 
dVes{4>,T)/dcj). In section HI, this derivative enters 
in the computation of the bounce solution of the eu¬ 
clidean equations of motion. For the numerical values 
of the SM input parameters we take Mw = 80.384 GeV, 
Mz = 91.1876 GeV, and v = 246.22 GeV (respectively, 
the W and Z pole mass and the vacuum expectation 
value of the Higgs field). For the pole Higgs mass we take 
Mh = 125.09±0.24 GeV according to the latest combina¬ 
tion of both ATLAS and CMS experiments [17]. Finally, 
for the strong coupling constant evaluated at Mz in the 
MS scheme (simply as hereafter) and the top quark pole 
mass we take as = 0.1184 ± 0.0007, Mt = 173.34 ± 0.8. 
As in [1], the latter is a naive combination of ATLAS, 
CMS and TeVatron measurements plagued by unavoid¬ 
able systematic error due to complicated Monte Carlo 
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modeling [18]. We will come back to this point in sec¬ 
tion IV. 

In fig. 1, notice that at T = 0 there is no electroweak 
minimum since we neglect the quadratic part in the tree 
level potential (see eq. (Al)). This approximation is well 
justified since we are interested in large field values. The 
potential at T = 0 exhibits the expected behavior, chang¬ 
ing sign around field values (j) ~ 10^°-10^^ GeV; this is 
the instability scale at which the quartic coupling A((^) 
crosses zero in its RG evolution, and the effective poten¬ 
tial develops the true vacuum. At T 7 ^ 0 (for definite¬ 
ness, we take in fig. 1 T = 10^® GeV) thermal corrections 
dominate over the T = 0 part until ^ « T; for </> ;g> T, 
on the contrary, they are exponentially suppressed (see 
eqs. (All,A12)), and therefore subdominant if compared 
with the T = 0 contributions. The shape of the effective 


M, = 173.34 GeV, = 125.09 GeV, ag = 0.1184 



0/T 


FIG. 2: SM effective potential (normalized with respect to 

four times its maximum value) at T = 10^® GeV. The field 
(j> scales as 4>/T. The values of the input SM parameters are 
shown in the plot label. In the insert, we zoom in the region 
close to the maximum (axis labels as for the outer plot). 


potential at finite temperature can be better visualized 
in fig. 2 where we show the effective potential, normal¬ 
ized with respect to its maximum value, as a function 
of the Higgs field rescaled according to the ratio ^/T. 
The effective potential changes sign at about 0 ~ ST; 
thenceforth, it sinks towards the true vacuum of the the¬ 
ory. Notice that the latter turns out to lie at extremely 
large field value, (p « 10^° GeV [19]. However, this is 
not a problem as soon as one assumes the SM to be valid 
up to the Planck scale: what really matters in terms of 
tunneling probability—at finite temperature as well as at 
T = 0—is the turning point of the bounce solution rather 
than the precise location of the true vacuum. The for¬ 
mer, as we shall clarify in the next section, never exceeds 


in our analysis Planck-scale values. 


III. BOUNCE SOLUTION AND THERMAL 
TUNNELING 


The vacuum decay in a scalar field theory with a po¬ 
tential characterized by an absolute minimum (the true 
vacuum) and a higher local minimum (the false vacuum) 
was first described in [2, 3]. The decay proceeds via a pro¬ 
cess called bubble nucleation, that is the tunneling from a 
false vacuum field configuration to a field configuration— 
the bounce—containing a bubble of approximate true 
vacuum embedded in a false vacuum background. 

At zero temperature the bounce <(>_b(?') is implicitly de¬ 
fined by the differential equation [2, 3] 


cPp 3 d(j) 


dr"^ 


r dr 


dVesjp) 

dp 


lim p(r) = 0 , 

r—xx) 


dp 

dr 


= 0 


r—0 


( 5 ) 

with r^ = -|- jr]^, r euclidean time. It corresponds 

to a field configuration that sits in the false vacuum at 
a long euclidean time ago (t — >■ — 00 ), and emerges at 
rest on the other side of the barrier at time t — 0 [2, 3]. 
The euclidean action for the 0(4) spherically symmetric 
solution of eq. (5) is 


5'E[(('(r)] = 27r^ 



1 

2 



+ Veffip) 


( 6 ) 


Let us now focus for the moment on the best-fit values 
Mh = 125.09 GeV, = 173.34 GeV, «« = 0.1184. In 
the left panel of fig. 3 we show the SM bounce solution 
obtained using the tree level RG improved quartic po¬ 
tential (red-dashed line) and the full two-loop effective 
potential (blue solid line). 

At finite temperature the bounce solution Pb{t') is im¬ 
plicitly defined by [4-6] 


dPp^2 dp 
dr"^ r dr 


dVes{P,T) 

dp 


lim p(r) = 0 , 

r—>-oo 


dp 

dr 


= 0 

r =0 


( 7 ) 

with now r = [f]. The euclidean action for the 0(3) 
spherically symmetric solution is 


‘S's [</'(?')] = dTT 



1 /#\ 
2\drJ 


Ves{P,T) 


( 8 ) 


In the right panel of fig. 3 we show the SM bounce solu¬ 
tion at finite temperature T = 10^^ GeV. Both at zero 
and finite temperature we solve numerically the bounce 
equation by means of the shooting method, without any 
kind of approximation for the effective potential.^ Two 
comments are in order. 


^ In appendix A we compare the full numerical result with the 
approximate solution at finite temperature proposed in [5], and 
often used in the literature. 
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M, = 173.34 GeV, M;. = 125.09 GeV, = 0.1184 



X = Y Mp 


M, = 173.34 GeV, M* = 125.09 GeV, as = 0.1184 






X = Y Mr 


FIG. 3: Bounce ai T = 0 (left panel) and T = 10^^ GeV (right panel). At zero temperature we show the bounce solution obtained 
considering a simple tree-level, RG-improved potential (dashed red line) and the full two-loop expression (solid blue line). The 
field (j) and the four-dimensional euclidean distance r are rescaled using the Planck mass Mp = 1.22 x 10^® GeV. At finite 
temperature we rescale the field as 4>(r) = x ip{r), with M^f, = 10 x T. We rescale the three-dimensional distance according 
to X = r X Mp, with Mp = V^max/0max. This prescription greatly improves the efficiency of the numerical shooting method 
used to solve eg. (7). In the insert, the red arrow pictorially indicates the bounce solution describing the thermal tunneling. The 
tip of the arrow corresponds to (^s(O) = 8.8 x 10^^ GeV. The values of the input SM parameters are shown in the plot label. 


At zero temperature and at the tree level, i.e. consid¬ 
ering the potential V{(p) = A0^/4, the bounce solution 
has the following analytical form 




[V R 

^\X\R^+r^ ’ 


(9) 


where R is an arbitrary scale reflecting the scale invari¬ 
ance of the potential. This degeneracy is broken by 
quantum corrections [20], and only one specific value of 
R —the one saturating the path integral, and dehning 
the size of the bounce, Rm in the following—is singled 
out. We can use eq. (9) to check the reliability of our 
numerical shooting method (dashed red line in the left 
panel of fig. 3). First, defining the size of the bounce via 
4>b{Rm) = <(>b(0)/2, we extract Rm = 434.33 x M(f^. 
Second, plugging back this number into eq. (9), and 
choosing pL = I/Rm for the renormalization scale in 
A(/r) [19, 20], we indeed find an exact match between 
our numerical solution and the actual bounce in eq. (9). 

At finite temperature we rescale—in order to improve 
the efficiency of the numerical shooting algorithm—the 
field as (j){r) = x with = 10 x T, and the 
three-dimensional distance according to a; = r x Mp, 
with Mp = x/foMAX/(('MAX ■ Throughout our analysis we 
always find the relation (/)b( 0)/T ^ 10 (see hg. 3, right 
panel, for the specific case with T = 10^^ GeV). This is 
the value of the field configuration at which the bubble 


of true vacuum is nucleated. The red arrow in the insert 
plot in the right panel of hg. 3 pictorially represents the 
bounce solution, with the tip at (/'_b(0). 

The vacuum decay rate per unit volume at hxed tem¬ 
perature T is [4-6] 

r(T) ~ I ^-S3[4>B{r)]/T ^ 


where Ep = S 3 [(j)B{r)] represents the energy of a bubble 
of critical size. In the left panel of hg. 4 we show the 
euclidean action of the bounce solution (j)B(Y) as a func¬ 
tion of the temperature for the best-ht values of Mu, Mt, 
and Og. The differential decay probability of nucleating 
a bubble at a given temperature T is given by [7] 


dP 

d\nT 



( 11 ) 


with To ~ 2.35 x 10“"* eV and ry the age of the Universe. 
Notice that this formula is valid only in a radiation- 
dominated Universe. In the right panel of hg. 4 we show 
the differential probability dP/dlogi^T as a function of 
the temperature. The total integrated probability is de- 
hned as 


rTut-o 


dP{r) 

dT' 


dT' . 


( 12 ) 


F’(Fcut-off) = 
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FIG. 4: Left panel. Euclidean action of the bounce solution (f>B{r) as a function of the temperature. Right panel. Plot of the 
differential probability dP / dlog^^T as a function of the temperature. In the insert, we show the value of as a function 

of the temperature. For a given cut-off scale (for instance, A = 10^® GeV, solid horizontal magenta line) the integration of 
dP lo^^Q P must be cut-offed at the temperature satisfying the condition ((b (ff) A (in this exainple T'cut —off 10^® GeV, 
vertical dashed magenta line). The values of the input SM parameters are shown in the plot label. 


Tcut-off is the cut-off temperature obtained imposing the 
condition = A, where A is the cut-off scale of the 

SM, for the moment assumed to be A = 10^® GeV. In 
the insert plot in the right panel of fig. 4 we show the 
values of 0 b(O) at different temperatures. The cut-off at 
A = 10^® GeV corresponds to a maximum cut-off value 
on the temperature Tcut-off — 10^® GeV, as expected 
since ^b(0)/T oo 10. Larger values of (fsiO) = A would 
correspond to a Planck-scale dominated tunneling tran¬ 
sition. The cut-off temperature plays a fundamental role 
in connection with the thermal history of the Universe. 
In section IV B we will discuss this aspect in detail. For 
the moment, in order to keep the discussion as simple 
as possible, we stick to the value A = 10^® GeV. In¬ 
tegrating the differential probability using eq. (12), we 
find P(rcut-off) = 5.22 X <C 1. Gonsequently, we 

conclude that the electroweak vacuum of the SM for the 
present central values of M^, Mt, and Og is unstable but 
sufficiently long-lived if compared to the age of the Uni¬ 
verse, even including thermal corrections with the highest 
cut-off scale A = 10^® GeV.® 


® Notice that the probability defined in eq. (12) is not normalized 
to one, and therefore—strictly speaking—it can not be inter¬ 
preted as a probability in the usual sense. The correct interpre¬ 
tation was given in [21] where it was shown that in the bubble 
nucleation process the fraction of space in the false metastable 
vacuum configuration is given by /false = while the frac¬ 

tion of space in the configuration of true vacuum is given by 


The total probability computed before turns out to be 
much larger than the corresponding probability evaluated 
at T = 0, that is ~ 10“®°° [19]. Said differently, the elec¬ 
troweak vacuum is still metastable but thermal correc¬ 
tions greatly enhance the tunneling probability. It sim¬ 
ply implies that—extending the previous computation to 
different values of Mh, Mt, and as in the allowed exper¬ 
imental ranges—the resulting instability bound will be 
much more stringent if compared with the one obtained 
at T = 0. A comprehensive analysis in the context of 
the phase diagram of the SM will be carried out in sec¬ 
tion IV. For the moment, as a warm-up discussion, let 
us now try to change only the value of Mt. In fig. 5 
we show how the total probability of thermal tunneling 
changes as a function of Mt for three different values 
of the Higgs mass, Mh = 124.0,125.09,127.0 GeV, with 
as = 0.1184. The total probability increases going to¬ 
wards larger values of Mt , and smaller values of Mh ■ For 
illustrative purposes, we show the region corresponding 
to the best-fit, 1- and 3 -ct confidence regions of Mt ac¬ 
cording to Mt = 173.34 ± 0.8 GeV. For Mh = 125.09 
GeV, we find that the total probability of thermal tun¬ 
neling equals one for values of Mt extremely close to the 
1-cr confidence region. This is a remarkable result, given 
that at r = 0 the instability bound is reached only for 


/true = 1 — e ^. Therefore if P <C 1 (F 1) all the space is in 
the false (true) vacuum field configuration. 
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as = 0.1184 



FIG. 5: Plot of the logj^Q of the total probability as a function 
of the top mass for three different values of the Higgs mass, 
Mh = 124.0,125.09,127.0 GeV. For illustrative purposes, the 
vertical blue lines mark the best-fit, 1- and 3-cr values accord¬ 
ing to Mt = 173.34 ± 0.8 GeV. 


verifies the condition P ^ 1, where the thermal 
tunneling probability is given in eq. (12). 

o At T = 0, the instability region (marked by the 
dashed red line) corresponds to a zero-temperature 
tunneling probability 


Vv 

p = exp 


r Stt^ 1 

. 3|A(m)L 


> 1 , 


(13) 


where ru is the age of the Universe and V\j t^. 

o In the metastability region (yellow) Aes (0) does be¬ 
come negative below the Planck scale, and the effec¬ 
tive potential develops a second minimum deeper 
than the electroweak one. However, the decay 
probability verifies P < 1. 


In section IV A— as a natural continuation of what al¬ 
ready discussed in section III —we show the phase dia¬ 
gram of the SM at finite temperature with the highest 
cut-off A = 10^® GeV, corresponding to Tcut-off — 10^® 
GeV. In section IV B we discuss the role of Tcut-off in the 
early Universe, thus assessing under which conditions the 
instability bound at finite temperature applies. 


A. Instability bound at finite temperature 


Mt > 178 GeV. Motivated by this result, we turn atten¬ 
tion to the full phase diagram of the SM. 


IV. THE PHASE DIAGRAM OF THE 
STANDARD MODEL AT FINITE 
TEMPERATURE 

The phase diagram of the SM is divided in three re¬ 
gions describing absolute stability, metastability, and in¬ 
stability of the electroweak vacuum depending on the val¬ 
ues of the SM parameters. Among them, the top mass, 
the Higgs mass, and the strong coupling at weak scale 
play a dominate role. At finite temperature, we add a 
fourth region in order to discriminate between instabil¬ 
ity at T = 0 and thermal instability. All in all, the four 
regions are defined as follows. 

o The absolute stability region (green) verifies the 
condition Aeff((/>) ^ 0 all the way up to the Planck 
scale.The effective potential does not develop a 
second, deeper minimum, and the electroweak vac¬ 
uum is stable. 

o The instability region at finite temperature (red) 


^ Aeff is the effective quartic coupling accounting for one- and two- 
loop corrections which is extracted from the RG-improved effec¬ 
tive potential. 


In fig. 6 we show the phase diagram of the SM in terms 
of the Higgs and top mass. The gray ellipses refer to the 
I-, 2-, and 3 -(t confidence regions obtained considering 
Mt = 173.3 ± 0.8 GeV and Mh = 125.09 ± 0.24 GeV. 
At T = 0, the instability bound correctly reproduce the 
known result [1] according to which, for instance, values 
^ 178 GeV are excluded if Mh — 125 GeV. 

At finite temperature, the situation drastically 
changes. As expected, the instability bound is pushed 
towards lower values of Mt- To fix the ideas, values 
^ 174.5 GeV are excluded if My — 125 GeV. Includ¬ 
ing the uncertainties on the strong coupling at the weak 
scale (dot-dashed lines in fig. 6) the bound becomes even 
more stringent, and values Mt > 173.6 GeV are excluded 
if Mh ~ 125 GeV and = 0.1163. 

At finite temperature, and assuming the highest cut-off 
scale A = 10^® GeV, the instability bound excludes, tak¬ 
ing into account the present experimental uncertainties 
on as , almost one half of the allowed experimental range 
for {Mh,Mt). In terms of Mt we extract the following 
bound 


< 174.459 -h 0.4285 x 

GeV 



-f 0.283 X 


-0.1184\ 
0.0007 J 


125.09 


(14) 


In fig. 7 we show the phase diagram of the SM in terms of 
the top mass and the strong coupling at the weak scale, 
keeping Mh fixed at Mh = 125.09 GeV. As before, the 
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FIG. 6: SM phase diagram at finite temperature and cut-off 
scale A = 10^® GeV. Solid (dashed) red line: instability bound 
with (without) thermal corrections. We also show the 1-, 2- 
, and 3-a ellipses corresponding to Mt = 173.3 ± 0.8 GeV 
and Mh = 125.09 ± 0.24 GeV (assuming a two-dimensional 
gaussian distribution without correlations). 


Ml, = 125.09 GeV 



Mt [GeV] 

FIG. 7: Same as in fig. 6 but in the plane {Mt,as). We 

also show the effect of a 1 GeV shift in the determination of 
the top pole mass (dashed ellipses). The Hiqqs mass is fixed 
at Mh = 125.09 GeV. 


ellipses mark the 1-, 2- and 3 -(t confidence regions with 
Mt as in fig. 6 and Us = 0.1184 ± 0.0007. For illustrative 
purposes, we also show (dashed ellipses) the effect of a 
1 GeV shift in the determination of the top pole mass. 


Such shift symbolically represents the systematic error 
involved in the naive combination of ATLAS, CMS and 
TeVatron results used in this paper, Mt = 173.34 ± 0.8 
GeV. Moreover, one should always keep in mind that 
the experimentally measured top mass is not the pole 
mass entering in the computation of the instability bound 
but the outcome of a complicated reconstruction of top 
quark decays (often dubbed the Monte Carlo mass). This 
fact amounts to a further source of uncertainty. As well 
known, and emphasized in this plot, the measurement of 
the top quark pole mass plays a crucial role in the deter¬ 
mination of the actual position of the SM in the phase 
diagram [22]. With the inclusion of thermal corrections, 
the situation becomes even more severe if compared with 
the r = 0 case, since now a small shift of the measured 
values can drastically change the phase of the electroweak 
vacuum in both directions, towards the stability as well as 
the instability region. This result motivates the need of 
a future high-energy electron-positron collider, where Mt 
could be unambiguously measured with a few hundred 
MeV accuracy in the scattering process e+e“ —>■ tt [18]. 

B. Instability bound and reheating temperature 

Thermal corrections are computed assuming the Higgs 
field in equilibrium with a thermal bath at temperature 
T. The occurrence of this condition strongly depends on 
the thermal history of the Universe. During inflation [23] 
all the energy is stored in the inflaton field, which slowly 
rolls down towards the minimum of its effective potential. 
Once reached, inflation ends, and the inflaton begins to 
oscillate near the minimum. SM particles are created be¬ 
cause of their interactions with the inflaton field: the ki¬ 
netic energy of the oscillating inflaton is gradually trans¬ 
ferred into the ultra-relativistic SM particles produced 
in the final state of its decay. Eventually, SM particles 
reach a state of thermal equilibrium at the temperature 
Tr,h, dubbed reheating temperature [24]. Thenceforth, 
the temperature scales according to T oc a“^, as in the 
ordinary radiation-dominated phase (as customary, a is 
the Friedmann-Robertson-Walker scale factor). Strictly 
speaking, the applicability of our computation is limited 
to T < Trh. In order to further investigate this impor¬ 
tant point, we explore two possible scenarios. 


1. Instantaneous reheating 

We start describing the reheating as an instantaneous 
process. In this case the decay probability is given by 
eq. (12), with Tcut-off = Trh. In fig. 8 we show how 
the instability bound changes for different values of Trh- 
As clear from the right panel of fig. 4, the largest con¬ 
tribution to the total probability comes from the high- 
temperature region, and a decrease in the cut-off quickly 
weakens the instability bound. We show the impact of 
different reheating temperatures in fig. 8. At Trh — 10^^ 









tally allowed. Moreover, the hypothesis of instantaneous 
reheating is a crude, yet not unrealistic, approximation. 
More likely, reheating is a dynamical process. In the next 
section we will elaborate this point and its consequences 
in more detail. 


2. Including the dynamics of reheating 

Reheating is not an instantaneous process. On the con¬ 
trary, the radiation-dominated phase at T < Trh follows 
a stage of matter domination during which the energy 
density of the Universe is dominated by the oscillations 
of the inflaton field [27, 28]. Temperature scales accord¬ 
ing to T oc [27, 28]; in other words, during the 

oscillating phase the Universe cools down more slowly— 
if compared with the scaling T oc of the radiation- 
dominated phase—because of the heating effect of the 
inflaton decay. As shown in [27, 28] the maximum value 
of the temperature is 


FIG. 8: Dependence of the instability bound on the reheat¬ 
ing temperature Trh, assuming instantaneous reheating. The 
orange dot-dashed lines correspond to different values Trh = 
10^ GeV, with—from top to bottom—x = 10,12,14,16. 


GeV the instability bound is pushed towards the border 
of the 3-cr band on {Mh,Mt). For smaller values of the 
reheating temperature, e.g. Trh = 10^° GeV, the SM 
reenters in the metastability region. The bounds in fig. 6 
are well described by the following parametric formula 


Mt 

GeV 


< 0.283 X 


/as - 0.1184\ 
0.0007 ) 


-f Cl X 


Mh Trh 

^ « X log.„ — 


_C3_ 

C4 X logio + C5 


(15) 

with Cl = 0.4612, C 2 = 1.907, C 3 = —1.2 x 10^, C 4 = 
—0.323, C 5 = —8.738. In concrete, taking My = 125.09 
GeV, as = 0.1163 (close to the 3-cr lower bound), and 
Trh = lOi® GeV we find Mt < 173.65 GeV. 

Before proceeding, let us pause for a moment to com¬ 
ment about the current experimental limits on the re¬ 
heating temperature. Despite its relevance in our under¬ 
standing of the early Universe, very little is known about 
the actual value of the reheating temperature. An obvi¬ 
ous lower bound can be obtained requiring a successful 
Big Bang Nucleosynthesis, and it turns out to be Trh > 
10 MeV [25]. As far as the upper bound is concerned, 
it is possible—assuming instantaneous reheating—to re¬ 
late the reheating temperature to the energy scale of the 
inflationary potential [26]; since the latter can be con¬ 
strained using the limit on the tensor-to-scalar ratio of 
the amplitudes produced during inflation, it is possible to 
extract a bound on Trh. All in all, one finds Trh < 10^® 
GeV [26] . High values of reheating temperature—as large 
as the ones considered in fig. 8 —are therefore experimen- 


Tmax = 



1/8 


gy^(TRH) 

^^"^(Tmax) 




(16) 


where g* (T) is the effective number of degrees of freedom, 
and Hf is the Hubble parameter at the end of inflation. 
The situation is schematically summarized in fig. 9. In 


1/^MAX T oc 1/^RH T oc a~^ IfT 

—I-1-► 

Oscillating phase Radiation—dominated phase 


FIG. 9: Schematic representation of the thermal evolution of 
the Universe after inflation. At the end of the reheating pro¬ 
cess (T < Tycr) the temperature scales according to T oc a~^, 
as in the ordinary radiation-dominated phase. During the os¬ 
cillating phase of the inflaton, before reheating is completed, 

roca-3/«. 


the region Trh ^ T ^ Tmax we can not compute the 
decay probability using eq. ( 11 ), since it relies on the 
assumption of a radiation-dominated Universe. Using 
the scaling T oc in the region Trh ^ T < Tmax 

the differential decay probability becomes [7] 




(17) 


All in all, the total integrated probability is given by 

dP{T') 


f 

T(TRH,i7/) = / 
Jo 


dT 


dT' (18) 


eq. (11) 


r^MAX ^p{T') 


dT' 


dT' 


eq. (17) 


and it depends on the reheating temperature and the 
value of the Hubble parameter at the end of inflation 
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via eq. (16). Notice that, for a given Trh, the Hubble 
parameter is characterized by the lower bound i/™™ = 
[ 47 r^( 7 *(TRH)/ 45 ]^/^(T^jj/Mp); this bound follows from 
the limit in which the inflaton energy density equals the 
energy density of a thermal bath with temperature Trh- 
In fig. 10 we show how the instability bound changes for 
different values of Trh including the dynamics of reheat¬ 
ing. For definiteness, we take Hf = 10^"^ GeV. As ex- 



FIG. 10: Dependence of the instability bound on the re¬ 

heating temperature Trh, including the dynamics of reheating. 
The magenta dot-dashed lines correspond to different values 
Trh = 10^ GeV, with—from top to bottom—x — 10,12,14. 
We take Hf = lO^'* GeV and g^T) = 106.75. 

pected, comparing the same values of the reheating tem¬ 
perature analyzed in fig. 8, the instability bound becomes 
more stringent including the dynamics of reheating. As a 
benchmark example, the value Trh = 10^° GeV—outside 
the experimental ellipses in fig. 8 — approaches again the 
edge of the 3-a region if the oscillating phase is included. 
In order to better investigate the role of the interplay 
between the reheating temperature and the Hubble pa¬ 
rameter at the end of inflation, in fig. 11 we recast the 
instability bound in the plane {Hf, Trh) for different val¬ 
ues of the top mass. For each value of Mt, the values of 
Trh above the corresponding red curve are excluded. We 
notice that the instability bound, for a fixed value of Mt, 
becomes stronger increasing the value of Hf; this is ex¬ 
pected, since the larger Hf the higher Tmax- However, 
we also notice that the Hf dependence is very mild (af¬ 
ter all Hf enters only as in Tmax)- As for the rest, 
fig. 11 retraces what already foreseen in fig. 10. Stringent 
bounds on the top mass—close to the present experimen¬ 
tally measured central value—can be reached only for 
very high (yet reasonable) reheating temperatures. For 
reheating temperatures Trh — 10^°-10^^ GeV, the bound 


M,, = 125.09 GeV, = 0.1184 



FIG. 11: Instability bound in the plane (H/jTrr), for dif¬ 
ferent values of the top mass. We include the dynamics of re¬ 
heating, and we keep fixed Mh = 125.09 GeV and as = 0.1184. 
The gray region is excluded by the condition Hf < For 

each Mt, the region above the corresponding red curve is ex¬ 
cluded. 


on the top mass is Mt > 176 GeV, at the border of the 
experimental S-cr confidence interval. 

Let us now conclude this section summarizing in a nut¬ 
shell our results. Thermal corrections are relevant for the 
computation of the instability region in the SM phase di¬ 
agram, and they can put a very stringent bound on Mt 
close to the present measured central value if also the un¬ 
certainties on Os are included. However, they crucially 
depend on the temperature of the early Universe. As al¬ 
ready noticed in [7, 8 ], therefore, the fate of the SM and 
its cosmological history are inextricably linked. 

A crucial question now seems to be: what was the high¬ 
est temperature ever recorded in the early Universe after 
inflation ended? On a general ground, one could be in¬ 
clined to think that it must have been very high. Let us 
provide one example in the context of thermal leptogen- 
esis and neutrino mass generation via type-I seesaw [7]. 
On the one hand, in order for baryogengesis to proceed 
via leptogenesis the mass scale M of the sterile neutrinos 
must be of the order of 10® GeV or larger [29, 30]; on 
the other one, in order to produce thermally the heavy 
neutrino states a reheating temperature of the Universe 
after inflation of Trh > M is required. This simple ar¬ 
gument seems to point towards a value of the order of 
Trh 10^® GeV, a temperature high enough to gener¬ 
ate large thermal corrections, as shown in fig. 10 . 

Moreover, as already stated in the introduction, in [ 8 ] a 
large reheating temperature after inflation (from Trh — 
10^ GeV up to Trh — 10^^ GeV, the actual value de- 
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pending on the instability scale of the Higgs potential 
and the value of the Hubble constant during inflation) 
seems to be suggested by inflation itself, since it may 
tame dangerous quantum fluctuations of the Higgs field. 


V. CONCLUSIONS AND PROSPECTS 

In this paper we revisited and updated the computa¬ 
tion of thermal corrections to the stability of the elec- 
troweak vacuum in the SM. We followed the approach 
of [6], based on i) the computation of the effective po¬ 
tential at finite temperature, and ii) the exact numeri¬ 
cal solution of the bounce equation. Although the im¬ 
portance of thermal corrections was recently reiterated 
in [7, 8, 31], a full computation including the most up¬ 
dated expressions for effective potential, beta functions 
and matching conditions was still missing. Our results 
can be summarized as follows. 

First, we studied the impact of thermal corrections 
on the instability of the electroweak vacuum consid¬ 
ering the highest allowed cut-off for the temperature, 
Tcut-off — 10^® GeV. The corresponding SM phase di¬ 
agram is shown in fig. 6 (referred to the parameters 
and Mt), and fig. 7 (referred to the parameters Mt and 
tts). Thermal corrections turn out to be very important, 
and they strengthen the constraining power of the insta¬ 
bility bound on the SM parameters if compared with the 
case at T = 0. If taken at face value, our results show 
that the instability bound at finite temperature excludes 
values of the top mass Mt > 173.6 GeV, if Mh — 125 
GeV, and including the uncertainties on the strong cou¬ 
pling constant at the weak scale. Parametrically, our 
bound is given by eq. (14). 

Second, we studied the temperature dependence of the 
instability bound. Thermal corrections crucially depend 
on the reheating temperature, hence on the cosmolog¬ 
ical history of the early Universe after inflation ended. 
From this perspective, the case previously studied cor¬ 
responds to a limit scenario in which Trh — 10^® GeV. 
In order to explore the temperature dependence, we in¬ 
vestigated two possible situations. 1) We considered the 
reheating after inflation as an instantaneous process. Ac¬ 
cording to this simplified assumption, the Universe expe¬ 
rienced a sharp transition from the inflationary epoch to 
the radiation-dominated phase. Our results are shown 
in fig. 8. The instability bound at finite temperature, 
now cut-offed at Tcut-off = TrHj weakens. However, for 
Trh — 10^^ GeV the instability bound still lies at the 
edge of the 3-cr confidence region for the experimentally 
measured values of Mh and Mt- For larger values of 
Trh, the SM enters in the instability region. Parametri¬ 
cally, our bound as a function of Trh is given by eq. (15). 
2) We included in our analysis the dynamics of reheat¬ 
ing. The instability bound becomes stronger if compared 
with the case of instantaneous reheating since it includes 
the oscillating phase of the inflaton field in the interval 
Trh T < Tmax, where Umax is given by eq. (16) and 


depends on the value of the Hubble parameter at the end 
of inflation. Our results are shown in fig. 10. We find that 
if Trh 10^® GeV the SM starts to fall in the instability 
region of the phase diagram. 

To conclude, the metastability region of the SM phase 
diagram considerably shrinks if thermal corrections to 
the decay of the electroweak vacuum are included. On 
the quantitative level, the impact of these corrections de¬ 
pends on the cosmological history of the early Universe, 
as shown in [7, 8] and discussed in more detail in this pa¬ 
per. From a more qualitative perspective, unveiling the 
true nature of near-criticality becomes an even more ur¬ 
gent question. To this end, possible directions include a 
better measurement of the top quark pole mass—if pos¬ 
sible at a future high-energy electron-positron collider— 
and a deeper understanding of the interplay with the 
physics of the early Universe. 
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Appendix A: Effective potential 


The effective potential is given by two contributions, 
the T = 0 corrections and the thermal effects, computed 
in the MS scheme and in the Landau gauge. For the zero- 
temperature term we have considered up to the two-loop 
corrections, but the complete expression is too lengthy to 
be given here. In order to setup our conventions, we only 
show the improved tree-level expression and the one-loop 
terms 


f4(<(') = + ~ ^A(t)(()^(t) (Al) 


ki_ioop(<(') — 'y ^ 


rij 




i=W,Z,t,x,h 

where the coefficients rii, Ci are 


rr^ ■ 


(A2) 


nw =6, nz = 3, n* = -12, = 3, Uh = 1, 

Cw=Cz = 5/6, Ct = C^ = Ch = 3/2, (A3) 
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while the mass parameters are given by 

mwit) = , (A4) 

^zit) = ^ : (A5) 

, (A6) 

mlit) = -^ « X{t^^t) , (A7) 

’mlii) = +3A(t)0^(t) « 3A(t)</)2(t) .(A8) 


Since we are interested in large field values, we neglect 
the quadratic term in the Higgs potential. All the SM pa¬ 
rameters are running with the three-loop RG equations, 
so that our analysis takes into account all the NNLL con¬ 
tributions. 

The one-loop thermal corrections to the effective poten¬ 
tial are (see [32] for a thorough discussion) 


Vi-ioopi4>^T) = 
Vnngi4>,T) = 


i^W,Z,x,h 


Jb 


rp2 


UtT^ (mlit) 
-Jf 


E 


i=WL,ZL,'lL,X,h 


127r 


mlit) 


rp2 


27r^ 

3/2 


2^2 


2^2 


3/2' 


(A9) 

(AlO) 


The thermal integrals are 


Jsiy) 

Jpiv) 


/ dxx^ In 

1 _ 

Jo 


/ dxx^ In 

1 + 

Jo 

L -1 


(All) 

(A12) 


The one-loop thermal potential is improved by the one- 
loop ring resummation of daisy diagrams in which only 
the bosonic degrees of freedom are taken into account 
and, in particular, only the longitudinal component of 
the vector fields. The degeneracy coefficients are 


nwL = 2, nz^ =1, = I . (A13) 

The Debye masses are Mli(p) = mlit) + Ilii(j),T), with 
the following temperature-dependent self-energies 

n.(<^,T) = + ^ + = ^ 

IlwA<l^,T) = 

IIwt{ 4',T) = UzT,i(j),T) =Il^^i(j),T) = 0 , (A14) 


where we omit the t-dependence implied by the RG im¬ 
provement. Finally, mapping (W3, B) into (Z, 7 ), we find 




mlit) 

mlit) 


11 

6 cos^ 
11 

6 cos^ 


T^ + Aicl),T) 
- Ai(l),T) 


(A15) 


with 


A 2 (^,T) 


:i't) + 


11 cos^ 29w 
3 cos^ 9^/ 


mlit) + 


11 
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12 cos^ 9w 


rp2 


X 

. (A16) 


Having set the formalism, let us now quantify the impact 
of the ring corrections in eq. (AlO). These corrections 
take into account the resummation of daisy diagrams. 
As clear from eq. (AlO), the ring contribution vanishes 
in the limit cj) ^ T. The numerical approach carried out 
in section HI showed that, for a given T, the thermal tun¬ 
neling always occurs at field value ^(> 5 ( 0 ) ^ 10 x T. This 
is enough to argue that ring contributions do not play a 
crucial role. Fig. 12— where we compute the total prob¬ 
ability in eq. ( 12 ) with (red) and without (magenta) ring 
contributions—confirms this hypothesis. Ring contribu¬ 
tion generate a ^ 5% correction to the bounce action 
(insert plot in fig. 12 ); in turn, this correction translates 
into a ~ 0.2 GeV strengthening of the instability bound. 

The analysis presented in this work is performed by 
numerical methods and, as such, does not rely on any 
analytical approximation. It is interesting, therefore, to 
compare our numerical results to those obtained, for in¬ 
stance, in the large T regime [5] . In the high-temperature 
limit the effective potential can be written in the form 


Vf^sicj), T) ~ + const (AI7) 
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Mh = 125.09 GeV, as = 0.1184 



M, [GeV] 


FIG. 12: Same as in fig. 6 but considering only Mh = 

125.09 GeV, and comparing the case with (red, filled circles) 
and without (magenta, filled sguares) ring corrections in the 
effective potential. In the insert, we compare the bounce action 
(same as m fig. f, left panel) with (solid, red) and without 
(dashed, magenta) ring corrections. 



My [GeV] 


FIG. 13: Same as fig. 6. The dot-dashed line corresponds to 
the instability bound in the high-temperature approximation. 

Appendix B: Beyond leading order thermal 
corrections 


where the constant term, being 0-independent, can be 
neglected and the coefficient is 

- 9 '^ + ^ 9 ^ + 8A -b . (A18) 

The first line comes from the high-T expansion of the 
thermal integrals while the last two from the ring po¬ 
tential. For large field values, the effective potential at 
T = 0 can be expressed in terms of an effective quartic 
coupling Aeff which accounts for one- and two-loop cor¬ 
rections. With such a simple expression for the effective 
potential, the bounce equation can be solved straightfor¬ 
wardly, obtaining [5] S'3[0B(r)] ~ — (6.015)7rK/AeffT. We 
recall once again that all the parameters are scale de¬ 
pendent and run with the RG equations. To minimize 
the impact of large logs, the renormalization scale is cho¬ 
sen to be equal to the canonical normalized scalar field. 
Moreover, as we have already shown, the thermal tunnel¬ 
ing is characterized by a field value roughly of the order 
of the temperature, thus reducing the analysis to a prob¬ 
lem with just one scale, fixed by the temperature. We 
show in fig. 13 the SM phase diagram within this high- 
temperature approximation. This results in a less tight 
instability bound of ^ 0.6 GeV with respect to the full 
numerical analysis. 


In this paper we truncated the perturbative expansion 
of the effective potential at finite temperature at one-loop 
(including resummed ring diagrams). In this appendix 
we discuss, at the qualitative level, the validity of this 
description together with possible future improvements. 
The effective potential at finite temperature enters in the 
euclidean action in eq. (8) . Since at zero temperature we 
worked at two-loop order, it is natural to ask what is the 
impact of two-loop thermal corrections. At two loops, 
thermal corrections to the effective potential were stud¬ 
ied in [33-35] in the context of the electroweak phase 
transition. Here what we want to stress is that adding 
two-loop thermal corrections to the effective potential at 
finite temperature does not improve the precision of the 
computation, since the one-loop result is already plagued 
by theoretical uncertainties—very likely of the same or¬ 
der of the two-loop corrections. The crucial point is that 
the euclidean action in eq. (8) relies on different approxi¬ 
mations. Before proceeding, we stress that a comprehen¬ 
sive analysis of the theoretical errors associated with the 
computation of the stability of the electroweak vacuum 
at finite temperature is an extremely difficult task—well 
beyond the purposes of this paper and, to the best of 
our knowledge, never studied before in the literature. In 
what follows, we highlight the most relevant aspects of 
such analysis. 

o High-temperature approximation. At sufficient 
large temperature, in the computation of the eu¬ 
clidean action the integration over the euclidean 
time amounts to multiply by T~^ the three dimen- 
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M,= 173.34 GeV,M;. = 125.09 GeV, = 0.1184 M, = 173.4 GeV,M;. = 125.09 GeV, = 0.1184 




FIG. 14: Euclidean action of the bounce solution as a function of the temperature. Left panel: we compare the case without (blue 
solid line) and with (red dashed line) the one-loop correction to the kinetic term, eq. (B2). Right panel: to test gauge-dependence, 
we compare the impact of different choices of ^ for the one-loop thermal corrections to the effective potential. 


sional action corresponding to the 0(3) symmetric 
bubble [36] 

SEl<(>B(r)]=T-^S3l<fB(r)] , (Bl) 

with S'3[0 b(7 ’)] as in eq. (8). The parameter con¬ 
trolling this approximation is the inverse of the 
bounce size at zero temperature [36] . In section III 
we computed this quantity, and we found ~ 
2.8 X10^® GeV. From the right panel in fig. 4, we see 
that the decay probability is dominated by larger 
values of temperature. However, the validity of the 
approximation in eq (Bl) is not always guaranteed 
and deserves further studies. 


o Corrections to the kinetic term. In the computation 
of the euclidean action one should also include— 
in addition to the one-loop thermal corrections to 
the effective potential—one-loop corrections to the 
kinetic term. In full generality, these corrections 
can be written in a gradient expansion in power of 
derivative of the classical background field cj) [37] 

[Hr, i)] = ^ ^ (OH) + Ves (H T) 
f °° 1 

+ / ’ (B2) 

/ ^ Tii 

"'P n=2 

where each Zn(4>, T), in turn, can be expanded in a 
power series in the couplings. The euclidean space- 
time integration is 


= / dt 




= l/T , 


(B3) 


o Gauge dependence. If the corrections to the ki¬ 
netic term are neglected, the effective action be¬ 
comes gauge-dependent as a consequence of a bro¬ 
ken Nielsen identity [38] . 

In fig. 14 we estimated the impact of gauge dependence 
and corrections to the kinetic term. In the left pane, we 
included the corrections to the kinetic term truncating 
at first order the gradient expansion in eq. (B2), and 
computing at one-loop the wave-function renormalization 
Z2(4>,T) following [37]. In the right panel, we estimated 
the impact of gauge dependence at finite temperature at 
one-loop in a generic gauge. In both cases we found 
a correction to the effective action of the bounce of the 
order of few percent. 

A more detailed analysis of these corrections will be 
presented in a forthcoming work. 
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